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Abstract 



We analyze the annual-modulation effect, measured by the DAM A Col- 
laboration with the new implementation of a further two-years running, in 
the context of a possible interpretation in terms of relic neutralinos. We im- 
pose over the set of supersymmetric configurations, selected by the annual- 
modulation data, the constraints derived from WIMP indirect measurements, 
and discuss the features of the ensuing relic neutralinos. We critically discuss 
the sources of the main theoretical uncertainties in the analysis of event rates 
for direct and indirect WIMP searches. 
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I. INTRODUCTION 



The effect of annual modulation measured by the DAMA Collaboration in its WIMP 
direct search experiment with a Nal(Tl) detector and reported in Ref. was analyzed in 
terms of relic neutralinos in Refs. [|]-^. In these papers, we proved that this interpretation 
is compatible with the DAMA data, and entails a relic neutralino which might have the 
role of a major component of dark matter in the Universe, especially when the uncertainties 
affecting the evaluation of the neutralino-nucleon cross section are taken into account 0]). 
We have also presented in detail other physical properties of such a neutralino, both in a 
Minimal Supersymmetric extension of the Standard Model and in supergravity schemes 
0, and we have outlined how indirect measurements of WIMPS (low-energy antiprotons in 
cosmic rays and up-going muon fluxes from the center of the Earth and from the Sun) may 
bring further information [0], by way of constraints on the supersymmetric configurations 
derived from the DAMA annual-modulation results 

New data, collected by the DAMA Collaboration in a further two-year running of the 
Nal(Tl) experiment for an exposure of 38 475 kg ■ day, and now presented in Ref. 0, confirm 
their previous finding of an annual-modulation effect, which does not appear to be related 
to any possible source of random systematics. Taking together all (old and new) samples of 
data for a total exposure of 57 986 kg ■ day, the effect turns out to be at a Aa C.L. Performing 
a maximum likelihood analysis in terms of and ^ o'f caiar ' where is the WIMP mass, 
'^scaiar'"'^'' the WIMP-uucleou scalar elastic cross-section, and ^ = Px/ Pi is the WIMP 
fractional amount of local non-baryonic dark-matter density pi, the DAMA Collaboration 
presents a 3cr C.L. annual-modulation region, in the plane "^x^^ '-'^scalar°°^ whose actual 
size depends on whether or not the upper-bound constraints previously obtained by the 
same Collaboration are included, and on the values assigned to the galactic astrophysical 
velocities. For the purpose of the analysis carried out in the present paper, among the regions 
presented in Ref. we select the one, which is obtained from the annual-modulation data, 
by including the upper-bound constraints of Ref. [0], by setting pi at the standard reference 
value: pi = 0.3 GeV cm~^, and by taking into account uncertainties in the astrophysical 
velocities of the usual galactic Maxwellian distribution (170 km s~^ < "Wq < 270 km s~^; Vesc 
= 450-650 km s^^; where vq is the rotational velocity of the local system at the position 
of the solar system and Vesc is the galactic escape velocity). This region is the one shown 
in Fig. 1 (should one include also a bulk rotation of the dark matter halo this region 
would elongate along the horizontal axis up to ~ 230 GeV @]). In this figure we also 
show the contour lines for the three values Vq = 170,220,270 km s~^, separately [^]. In 
the comparison of the experimental data with the theoretical evaluations one has to further 
consider the uncertainty in p/: 0.1 GeV cm~^ < p« < 0.7 GeV cm~^ |]10|,|Tl]]. Fig. 2 displays 
how the DAMA annual-modulation region shifts along the vertical axis, as the value of pi 
is varied within its uncertainty range. The four panels correspond to the representative 
values: pi = 0.1, 0.3, 0.5, 0.7 GeV cm~^. In Fig. 2, as well as in all subsequent figures, 
where experimental results of direct and indirect WIMP measurements are compared with 
theoretical evaluations, separate panels are used for the four representative values of pi. 

In the present paper we investigate the implications of the DAMA data with the total 
exposure of 57 986 kg ■ day in terms of relic neutralinos, along the lines previously developed 
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in Refs. [Q-Q. We single out the set of the supersymmetric configurations compatible with 
the DAMA data, and then apply to this set, denoted as set S, the constraints derived 
from experimental indirect searches for WIMPs (up-going muons at neutrino telescopes and 
antiprotons in cosmic rays). In this analysis we incorporate recent, and quite significant, 
theoretical and experimental developments. 

The supersymmetric theoretical framework adopted here is the Minimal Supersymmetric 
extension of the Standard Model (MSSM) [|T2|, which conveniently describes the supersym- 
metric phenomenology at the electroweak scale, without too strong theoretical assumptions. 
This model has been extensively used by a number of authors for evaluations of the neu- 
tralino relic abundance and detection rates (a list of references may be found, for instance, 
inH). 

The neutralino is defined as the lowest-mass linear superposition of photino (7), zino 
(Z) and the two higgsino states {H^, HI): x = ^il + ci2Z + a^H^ + aiH2. 

The MSSM contains three neutral Higgs fields: two of them {h, H) are scalar and one 
(A) is pseudoscalar. At the tree level the Higgs sector is specified by two independent 
parameters: the mass of one of the physical Higgs fields, which we choose to be the mass 
niA of the neutral pseudoscalar boson, and the ratio of the two vacuum expectation values, 
defined as tan/? = {H2)/{Hi). Once radiative corrections are introduced, the Higgs sector 
depends also on the squark masses through loop diagrams. The radiative corrections to the 



neutral and charged Higgs bosons, employed in the present paper, are taken from Refs. [|14 
The other parameters of the model are defined in the superpotential, which contains all 
the Yukawa interactions and the Higgs-mixing term ^HiH2, and in the soft-breaking La- 
grangian, which contains the trilinear and bilinear breaking parameters and the soft gaugino 
and scalar mass terms. 

To cast the MSSM, which originally contains a large number of parameters, into a form 
adequate for phenomenology, we follow the common procedure of introducing a set of re- 
strictive assumptions at the electroweak scale: a) all trilinear parameters are set to zero 
except those of the third family, which are unified to a common value A; b) all squarks and 
sleptons soft-mass parameters are taken as degenerate: mj, = rriq. = mo, c) the gaugino 
masses are assumed to unify at Mqut, and this implies that the U{1) and SU{2) gaugino 
masses are related at the electroweak scale by Mi = (5/3) tan^ 6]yM2. 

Once these conditions are implemented in the model, the supersymmetric parameter 
space consists of six independent parameters. We choose them to be: M2, /x, tan /?, uia, mo, A 
and vary these parameters in the following ranges: 10 GeV < M2 < 1 TeV, 10 GeV < < 
1 TeV, 80 GeV < < 1 TeV, 100 GeV < mo < 1 TeV, -3 < A < +3, 1 < tan/3 < 50. 
We remark that the values taken here as upper limits of the ranges for the dimensional 
parameters, M2, fi, mo, mA, are inspired by the upper bounds which may be derived for these 
quantities in SUGRA theories, when one requires that the electroweak symmetry breaking, 
radiatively induced by the soft supersymmetry breaking, does not occur with excessive fine 
tuning (see Ref. |jl5| and references quoted therein). 

We have further constrained our parameter space, by taking into account all the new 
experimental limits obtained from accelerators on supersymmetric and Higgs searches (LEP2 
p!6| , CDF [0). Notice that the new bounds from LEP2 and CDF constrain now rather 
severely the susy space, especially in the region of interest for direct detection (small ruh 
and, partially, large tan P 
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Moreover, the constraints due to the 6 — s + 7 process |]T8[ have been taken into account. 
In our analysis, the inclusive decay rate BR(i? — X^'y) is calculated with corrections up 
to the leading order. Next-to-leading order corrections |]19| are included only when they 
can be applied in a consistent way, i.e. both to standard-model and to susy diagrams. We 
require that our theoretical evaluation for BR(i? Xs'~f) is within the range: 1.96 xlO^^ < 
BR(i? —>■ Xs'y) < 4.32 xlO""^. This range is obtained by combining the experimental data 
of Refs. pPlpTj l at 95% CL. and by adding a theoretical uncertainty of 25%, whenever the 
still incomplete next-to-leading order susy corrections cannot be applied. 

Our parameter space has been further constrained by the request that the Lightest Su- 
persymmetric Particle (LSP) is the neutralino, rather than the gluino or squarks or sleptons. 
The current upper bound for cold dark matter may be establish as QcDAih'^ < 0.3 {h is the 
usual Hubble parameter, defined in terms of the present-day value Hq of the Hubble con- 
stant as h = ifo/(100 km s~^ Mpc~^)), on the basis of the most recent cosmological data 

. However, for sake of presentation of the results of the present analysis, which, anyway, 
never entail values of ^^^"^ excess of 0.6 (see last section), we do not impose the bound 
^CDuh^ < 0.3 in our selection of susy configurations. The neutrahno relic abundance is cal- 



culated here as illustrated in Ref. ||23|. We have checked that susy configurations which could 
potentially lead to coannihilation effects |Q are marginal in our selected super symmetric 
parameter space. 

A few comments are in order here. The restrictive assumptions a) - c) adopted above 
in the framework of the MSSM are instrumental in reducing the otherwise large number 
of independent parameters to a handful set of them (six in our scheme), and in making 
the calculations of a number of crucial observables (such as relic abundances and event 
rates) manageable. The few independent parameters of this simplified MSSM have the 
role of relevant scales for some fundamental quantities, such as scalar masses and gaugino 
masses, which in turn determine the size of the numerical outputs. This version of MSSM 
is obviously the simplest scheme for a susy model, and the most natural one to start with. 
However, one has to be aware of the fact that new experimental data could eventually force 
one to adopt more involved versions of supersymmetric models, for instance by relaxing 



some GUT-inspired relation (such as Mi ~ O.5M2) Ref. [^], or by including CP-violating 
phases |]26[| . 

As regards the distribution of relic neutralinos in our Galaxy, to start with we have as- 
sumed a standard halo population with a Maxwellian velocity distribution, whose dispersion 
speed is centered around 270 km s~^ (i.e., Vq = 220 km s~^). However, in the implementation 
of constraints from up-going muons at neutrino telescopes, we have also considered recent 
theoretical developments which may have quite contrasting effects on the expected signals 



27| , |28[| . These different instances are examined in Sect. HI. 

Data on antiprotons in space, combined with recent evaluations of the secondary an- 
tiproton component in cosmic rays due to spallation processes, are employed in Sect. IV 
to put further constraints on the original set S of susy configurations, singled out by the 
DAMA data. 

We give the results of our combination of the annual-modulation data with indirect 
measurement constraints in Sect. V, where we also discuss the cosmological properties for 
our set of relic neutralinos and present our conclusions. 
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II. SET OF SUPERSYMMETRIC CONFIGURATIONS SINGLED OUT BY THE 

ANNUAL MODULATION DATA 



In our papers of Ref. we proved that the DAMA annual-modulation region of Ref. 
HI] is widely compatible with an interpretation in terms of relic neutralinos, by showing 
that a sizeable portion of that region is covered by supersymmetric configurations, sat- 
isfying all accelerator bounds. Now we show in Fig. 1 that the new, more constrained 
annual-modulation region of Ref. ||^ is still largely compatible with the relic neutralino 
interpretation, though the supersymmetric space is now more severely constrained by the 
current limits from accelerators [ P^p!7| . 



In deriving the scatter plot shown in Fig. 1 we have used the scan of the susy pa- 
rameter space defined in the previous section. The neutralino-nucleon cross section has 
been calculated with the formulae reported in Ref. 0. As discussed in Ref. [Q, this cross 
section suffers from significant uncertainties in the size of Higgs-quark-quark and squark- 
quark-neutralino couplings. In fact, these couplings depend on quark masses rriq and quark 
scalar densities in the nucleon (qq) , which are still rather poorly determined. To be specific, 
we refer to the following quantities: the fractional strange-quark content of the nucleon 
y = 2 < ss > / {< uu + dd >) , the quark mass ratio r = 2ms/ {^u + nT'd), and the products 
mq < qq >'s. In our analysis we have taken into account the uncertainties in these quan- 
tities. Thus, our scatter plots comprise representative points which have been derived by 
using both of the two following sets of values, cumulatively: 

Set 1: y = 0.33, r = 29, 

mi <Il> = 23 MeV, ttIs < ss > = 215 MeV, mh<hh> = 50 MeV. (1) 



Set 2: y = 0.50, r = 29, 

mi <ll> = 30 MeV, < ss > = 435 MeV, nih < hh > = 33 MeV. (2) 

In Eqs. (|T]- ^ / stands for light quarks, s is the strange quark and h = c,b,t denotes heavy 
quarks. For the light quarks, we have defined mi <ll > = \[mu < uu > +mii < dd >]. Set 
1 and set 2 bracket, at least partially, the present uncertainties. In Sect. V.B, in connection 
with neutralino cosmological properties we will also mention the consequences of using a 
more extreme set of values (set 3 of Ref. 0). For the derivation of the values of the various 
sets see Ref. It is worth noticing that a new derivation of the pion-nucleon sigma term, 
c^TT AT, points to rather high values: cx^ at = 73.5±9 MeV [|^. By itself, this new result would 
increase the value of the quantity < ss > given in Eq.(^) by ~ 30%. We recall that the 
quantity < ss > is crucial in establishing the size of o"^ca/ar°"^ 12^ 



As for the values to be assigned to the quantity ^ = p^/ pi we have adopted a standard 
rescaling recipe. For each point of the parameter space, we take into account the relevant 
value of the cosmological neutralino relic density. When flx^"^ larger than a minimal value 
(n/i^)min, compatible with observational data and with large-scale structure calculations, we 
simply put ^ = 1. When Q^h'^ turns out to be less than (f2/i^)min, and then the neutralino 
may only provide a fractional contribution to dark matter, we take ^ = Q^h'^ /{Qh'^)yam- The 
value to be assigned to {flh'^)min is somewhat arbitrary, in the range 0.01 < (f2/i^)min < 0.3. 
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We use here the value (r2/i^)min = 0.01, which is conservatively derived from the estimate 

^galactic ~ 0.03. 

As we mentioned above, Fig. 1 shows that the annual-modulation region (here depicted 
for pi = 0.3 GeV cm~^) is largely covered by the scatter plot. This turns out to be the 
case also for the other representative values of pi, as is shown in Fig. 2. In each panel of 
this figure we only display the portion of the susy scatter plot which is contained in each 
of the relevant experimental region. In going from the generic scanning used for Fig. 1 to 
the one employed for Fig. 2, although keeping the overall range of variation of the susy 
parameter space, we have optimized the numerical scanning in order to have a number of 
configurations, large enough for our subsequent analyses. The covering by the scatter plots 
of the annual-modulation regions pertaining to different values of pi is more extended for 
large values of pi than for the small ones, as expected from the features of the generic plot 
of Fig. 1. 

We define as set S of susy configurations the set comprised of the configurations whose 
representative points in the plane fn-x^'^icaiar"^^ inside the annual-modulation regions 
displayed in Fig. 2. Only configurations of set 5* are retained in the analyses presented 
hereafter. We remark that set S is the union of all the subsets of susy configurations which 
refer to each of the following representative values for pi and Vq: pi = 0.1,0.3,0.5,0.7 GeV 
cm~^, Vq = 170,220,270 km s~^, separately. At any stage, our results will be analysed and 
presented in our figures in terms of the chosen representative values of pi and Vq, separately. 

Another experiment of WIMP direct detection, the CDMS experiment [^], is now en- 
tering the DAMA sensitivity region. The current CDMS upper bounds (either with or with- 
out subtractions) concern the upper left corner of the annual-modulation regions, which is 
anyway poorly populated by susy configurations; thus, they are currently marginal in con- 
straining the set S. The upper left corner of the annual-modulation regions is also partially 
disallowed by independent constraints due to indirect measurements (see Sect. V). 

Now we turn to the constraints which may be applied to the set S, using data from 
WIMP indirect search experiments. We set the limits for exclusion at the same CL. to 
which the DAMA region is currently set, i.e. 99.7% CL. 



III. CONSTRAINTS FROM NEUTRALINO NEUTRALINO ANNIHILATION 

INSIDE EARTH AND SUN 



Indirect evidence for WIMPs in our halo may be obtained at neutrino telescopes by 
measurements of the up-going muons, which would be generated by neutrinos produced by 
pair annihilation of neutralinos captured and accumulated inside the Earth and the Sun 
3^^4[. The size of the expected muon fluxes strongly depends on how these relic particles 



are distributed in the phase space and on the intrinsic efficiency of the celestial body in 
capturing the surrounding WIMPs. 

In the case of the Sun the capture rate is essentially determined by its strong gravitational 
field and by the size of the cross section of neutralino scattering off single protons. Instead, 
in the case of the Earth the capture process may quite significantly be enhanced by coherent 
neutralino-nucleus cross sections, whose size depends on mass-matching condition between 
and the nuclear mass of the dominant chemical constituents of the Earth (O, Si, Mg, 
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Fe) H. 

As for the phase-space neutrahno distribution in our neighbourhood, together with the 
usual one based on the standard Maxwelhan velocity distribution, whose dispersion speed is 
centered around 270 km s~^, we also consider two intriguing and conflicting models which 
have been recently discussed in the literature. 

Damour and Krauss [^] have proposed the existence of a solar-bound population, with 
velocities restricted to rather low values, < 50 km s~^ (for other papers on hypothetical 
solar-bound WIMP populations, see Refs. p6|-^). The Damour-Krauss solar-bound pop- 
ulation would have been produced by WIMPs which scattered off the Sun surface and were 
set (by perturbations from other planets) into orbits which cross the Earth orbit, but not 
the Sun. The ensuing velocities would be distributed in the range 25 km s~^ < f < 50 km 
s~^. This population, although totally irrelevant for the direct measurements by the DAMA 
Nal-detector, whose electron-equivalent threshold energy is 2 keV, has been shown to be 
potentially important in making the capture of relic neutralinos by the Earth particularly 
efficient, with a consequent enhancement of the expected output of up-going muons from 
the Earth, as compared to the standard one . For simple kinematical reasons, the lower- 
speed cut off implies that this enhancement is limited to WIMPs of masses lower than ~ 
150 GeV. 



On the other side, Gould and Alam [2S], using arguments based on calculations of as- 



teroids trajectories |41], have pointed out that solar-bound WIMPs could evolve in a way 



quite different from the one derived in Ref. [27], with an ensuing suppression of the up- going 



muon flux usually expected from the center of the Earth for a standard halo population. 
This suppression would be significant for WIMP masses above ~ 65 GeV. 

In the present paper we take into considerations all of these possible instances. First we 
consider the standard situation of a Maxwellian velocity distribution over the whole speed 
range, then we proceed to a critical examination of the other two cases, in which the low- 
speed interval is either overpopulated or de-populated ||2^ , as compared to the standard 
one. 

The neutrino spectrum and the ensuing up-going muon fiux are calculated as ex- 



plained in Refs. p3| , |3^ . Their normalization is set by the annihilation rate Ta of the neu- 
tralinos inside the celestial body (Earth or Sun), and depends, in turn, on the capture rate 
C of the relic neutrahnos by the celestial body through the formula F^ = ^tanh^ (^^^^ p2| , 
where t is the age of the macroscopic body (t ~ 4.5 Gyr for Sun, Earth) and = (CC^)"^^^, 
Ca is the annihilation rate proportional to the neutralino-neutralino annihilation cross- 
section and C denotes the capture rate. In a given macroscopic body the equilibrium be- 
tween capture and annihilation [i.e. F^ ~ C/2 ) is established only when t > ta- Whereas, 
in the case of the Sun, the capture-annihilation equilibrium is usually reached, due to the 
much more efficient capture rate due to the stronger gravitational field, for the Earth, the 
equilibrium condition is not easily realized. 

For the case of the standard halo population with a Maxwellian velocity distribution, 
C and Ta are calculated as in Refs. fS^JS^, and the ensuing muon fiux is denoted by 
j^^Earth^std p^^^ Damour-Krauss population the quantities C and F^ are evaluated 
according to the formulae of Ref. (the relevant muon fiux is denoted by (^<|)Earth-jDK-j p^^, 
the model conjectured by Gould and Alam PB|, we have applied to the standard capture rate 
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a suppression factor, which we have re-calculated ab initio in the scheme denoted as ultra 
conservative in Ref. 



28 



to cover the whole range of masses involved in the present paper. 
For many susy configurations the suppression factor in the ensuing up-going muon fluxes 
from the center of the Earth is stronger than the reduction factor in the capture rate alone, 
due to the relation between Ta and C, previously mentioned. For these configurations a 
reduction in the capture rate induces in the muon flux an extra suppression due to a critical 
increase in the time required for reaching equilibrium. The muon flux calculated in the 
Gould-Alam model is denoted here as (<|)Earth^GA_ 

All our neutrino fluxes include neutrino oscillations and use the procedure outlined in 
Ref. |^3|. Here we assume z/^ —>■ Uj. oscillations, with values for the oscillation parameters 
which are taken from the best fit performed in Ref. over the whole set of experimental 
data on atmospheric neutrinos: Am^ = 3 ■ 10~^ eV^, sin^ = 1. 

Some of our results are presented in Figs. 3-6, where we report various muon fluxes (or 
ratios of them) versus m^, for the four representative values of pi. The solid lines, depicted 
in Fig. 3 and Fig. 6, denote the 99.7% C.L. upper bounds, (<|)Earthyim^ derived from the data 
of the MACRO experiment from the center of the Earth and from the Sun, respectively 
(for similar limits from the Baksan experiment see Ref. ||46|| ). 

The scatter plots of Fig. 3 display some expected characteristic features, such as the 
peak at ~ 50-60 GeV, due to the mass-matching between and mpe- We notice 
that a number of configurations induce up-going muon fluxes in excess of the experimental 
bounds. Figs. 4-5 show what would be the enhancement or the reduction effect in in 
the case of the Damour-Krauss population or in the Gould-Alam conjecture, respectively. 
The size of these effects agree with the evaluations in Refs. P0| , p8| . For the Damour-Krauss 
population, the enhancement effect for some susy configurations appears larger here than 
in Ref. this is due to configurations (not considered in where rescaling in is 
effective. In Fig. 6 we display the scatter plots for the up-going muon flux from the Sun, 
expected for the standard halo population. The current experimental bound P3| sets quite 
marginal constraints. 

In Sect. V we use the results of this section to constrain the susy configurations of 
set 5*. The question, as of which model for the low-speed WIMP population among the 
two extremes of Refs. is applicable, is still open. Thus, we implement here the 

experimental bounds on the standard flux of up-going muons; namely, we exclude those 
configurations, whose (<|)Ea,rth-jstd -g excess of the 99.7% C.L. upper bound derived from 
the MACRO data. 



IV. CONSTRAINTS FROM COSMIC RAY ANTIPROTONS 



The possibility that annihilation of relic particles in the galactic halo might distort the 
spectrum of cosmic-ray antiprotons at low-kinetic energies {Tp <1 GeV) has been consid- 
ered by many authors ||47|-^. Indeed, in this energy range, the production of secondary 
antiprotons by interactions of primary cosmic-ray protons with the interstellar hydrogen has 
a kinematical drop off [0], which primary p's, created by relic neutralinos of appropriate 
mass and composition, might fill in, at least partially. The effectiveness of this argument 
to disentangle ordinary spallation contribution from a possible exotic component due to 
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relic particles depends dramatically on how accurately the secondary spectrum is calculated 

11,1113. 



This point was addressed in Ref. In that paper we improved the evaluation of the 
energy losses undergone by secondary antiprotons during their diffusion inside the Galaxy, we 
noticed that the as-yet most recent experimental data (BESS95 |5^) were fitted reasonably 
well by the secondary spectrum alone, and we examined critically how much room was 
still available, in the low-energy spectrum, for a contribution from an exotic component. 
Now, new experimental data (BESS97 |Q) and improved evaluations of the secondary 
spectrum |]5D| , |5^ , |5B[] further constrain the room left for primary sources. These instances, 
instrumental in making the separation between primary and secondary antiprotons more 
difficult, nevertheless confer to the cosmic-ray antiproton measurements a potentially more 
important role in establishing stringent constraints for relic neutralinos of relatively low 
mass in our halo, once some of the sizeable, still persisting, uncertainties are reduced. 

In the present work we have evaluated the primary antiproton flux, expected from neu- 



tralino annihilation, as in Ref. |4^, restricting the supersymmetric configurations to those of 
set S. We refer to for all the details concerning the evaluation of the production of these 
primary antiprotons as well as for the properties related to their propagation in the halo 
and in the heliosphere. Here we only recall the features of the neutralino mass distribution 



function adopted in [ESI as well as here. This mass distribution function is taken spheroidal 



and parameterized as a function p^i^, z) of the radial distance r from the galactic center in 
the galactic plane and of the vertical distance z from the galactic plane in the form 

where a is the core radius of the halo, Tq is the distance of the Sun from the galactic center 
and / is a parameter which describes the flattening of the halo. Here we take the values: 
a = 3.5 kpc, Tq = 8 kpc. In the case of a spherical halo (/ = !), we use the value pi = 0.3 
GeV cm~'^. When / < 1 (oblate spheroidal distribution), pi is taken as |]57| , pT) 



For each value of pi and of the relevant value of /: p//(GeVcm~'^) = 0.1 (/ = 1), 0.3 
(/ = 1), 0.5 (/ = 0.50), 0.7 (/ = 0.33), we have evaluated the top-of-atmosphere (TOA) 
antiproton fluxes, as the sum of the secondary flux and of the primary flux due to neutralino 
annihilation for the various supersymmetric configurations of set S, pertaining to that spe- 



cific value of pi- The secondary flux has been taken from Ref. [g^. Re-acceleration effects in 
the cosmic rays propagation, which might also be relevant for the features of the secondary 
antiproton spectrum at low energies [^,^], are not included here. Solar modulation has 
been evaluated according to the procedure discussed in Ref. [^. We have compared our 
theoretical results with the combined experimental data of BESS95 and BESS97 0], over 
the whole experimental energy-range (0.18 GeV < Tp < 3.56 GeV), by a calculation. 

The results are reported in Fig. 7. In the evaluation of the x^y addition to the ex- 
perimental errors, we have also taken into account the theoretical uncertainties, estimated 
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according to the results in Refs. p9| , |50| , with their appropriate energy dependence. Orien- 
tatively, they are in the following ranges: ± (45-55)% for the primary fluxes, ±(60-75)% 
for the secondaries, depending on the energy bin. 

In the following, we adopt the selection criterion of excluding from set S the configuration 
whose reduced is above the value xl = '^■4:4:, which corresponds to a 99.7% C.L. for the 13 
d.o.f. of the BESS 95+97 data. From Fig. 7 we notice that, especially at large values of pi, 
this constraint disallows a number of susy configurations. The reason why the cosmic-ray 
antiprotons constraint is not more effective in constraining set S is to be attributed mainly 
to the current large uncertainties affecting the evaluation of antiproton propagation in the 
galactic halo and in the heliosphere. 



V. RESULTS AND CONCLUSIONS 

Now we apply the experimental bounds from indirect searches discussed in Sects. III-IV 
to constrain the supersymmetric configurations of set 5*. 



A. Combining direct and indirect measurements 

Fig. 8 displays the extent of the covering of the annual-modulation regions (one for 
each value of pi) by the susy configurations, when the MACRO upper bounds are applied to 
^^Earth-jstd^ A comparisou of this figure with Fig. 2 shows that the implementation of these 
limits somewhat de-populate the covering regions, with a marked effect for the value of the 
neutralino mass which matches the mass of Iron, as expected. Apart from this, the extent 
of the regions covered by the scatter plots does not significantly change. 

Fig. 9 depicts what would be the effect for a solar-bound WIMP population a la Damour- 
Krauss. Especially at low values of pi there would be some shrinking of the original regions 
of the scatter plots in their upper parts, but still the annual-modulation regions would be 
widely covered by physical susy configurations. At variance with this case, the Gould-Alam 
conjecture would relax the consequences of the constraints applied in obtaining the plots of 
Fig. 8. 

Now, we return to the case where the experimental bounds (<|)Earthyim applied on 
j^^Earth^std^ When, ou top of these constraints, we also implement the constraints due to 
cosmic-ray antiprotons, we obtain that the scatter plots of Fig. 8 become somewhat de- 
populated, but without any appreciable modification in the contours of the covering regions, 
except for a quite marginal downward shift in their upper-left parts. Therefore Fig. 8 
may be considered as the final situation of our analysis, once also the implementation of 
the antiprotons constraints has been applied. We denote as set T the subset of S which 
comprises the susy configurations not disallowed by bounds on the standard up-going muon 
fluxes and on cosmic-ray antiprotons. 

We have analyzed the main properties of the configurations of set T; some of them are 
displayed in Figs. 10-11. We recall that the scatter plots of these figures are derived, as 
all previous ones, by using for the hadronic quantities, discussed in Sect. II, set 1 and set 
2, cumulatively. In Fig. 10 we note that the configurations of set T cover only a specific 
region of the susy parameter space not yet disallowed by accelerator constraints. The shape 
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of the distribution of the representative points of T in the plot of Fig. 10 is simply explained 
by the fact that the values of the scalar neutralino-nucleon cross section at the level of the 
DAMA data require either a large tan/3 or a small m/j (or both of these two conditions). 
This constraint is stronger when the values of the hadronic quantities are restricted to set 1, 
alone. Fig. 11 displays a correlation among uia and uiq which is mainly due to the interplay 
of these two quantities in generating a light nih- Again, restricting the scatter plot to points 
belonging to set 1, this correlation becomes more pronounced. We recall that, at variance 
with constrained sugra-supersymmetric models, in the MSSM we are using here, and 
mo are treated as independent parameters. 



B. Cosmological properties 

We turn now to an analysis of the cosmological properties of relic neutralinos of the 
susy configurations of set T. The relevant plots f^^/i^ vs are displayed in Fig. 12. It 
is remarkable that the region of main cosmological interest: VL^h"^ > 0.03 turns out 



to be widely populated, with values of fl^h"^ which approach, and even exceed, what may 



be considered as the current upper bound for cold dark matter: VtcDuh'^ < 0.3 [^. This 
means that the DAMA annual-modulation data are compatible with a neutralino as a major 
component of dark matter. We stress that the scatter plot would even shift upward, should 
we use for the hadronic quantities discussed in Sect. II the following set: y = 0.50, r = 36, 
mi <ll > = 33 MeV, nis < ss > = 585 MeV, ruh < hh > = 21 MeV. This set of values, 
denoted as set 3 in Ref. is more extreme as compared to set 1 and set 2, but still 
compatible with the current uncertainties. 

Finally, we notice that a rather strong de-population in the plots of Fig. 12 is present 
around Q^h'^ ~ 0.01 and for large values of pi. This effect is induced by the cosmic-ray 
antiproton constraint, since the calculated p fluxes have their maximal values for ^2^^^ close 
to the value below which we apply the rescaling of the local density, i.e. (f2/i^)min = 0.01. 
This property is quite general in this class of calculations, and it was already commented 
upon, for instance, in Ref. 0]. 



C. Conclusions 

In the present paper we have examined the possibility that the annual-modulation effect, 
measured by the DAMA Collaboration at a 4a confidence level 0, may be interpreted in 
terms of relic neutralinos. We have examined this problem, by employing the Minimal 
Supersymmetric extension of the Standard Model, as a model which does not impose too 
strong theoretical prejudices on the phenomenological analysis. We have taken into account 
all experimental constraints, from accelerators and from WIMP indirect experiments. 

Let us now summarize our main conclusions: 

• The annual-modulation effect mentioned above turns out to be compatible with an 
interpretation in terms of relic neutralinos. 
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• The set of supersymmctric configurations selected by the annual-modulation data is 
only modestly reduced by current experimental data from WIMP indirect searches 
(up-going muons from the Earth and the Sun, and cosmic-ray antiprotons) . 

• The set of supersymmctric configurations, selected by the annual-modulation data 
and not disallowed by the indirect measurements, comprise configurations of relevant 
cosmological interest, with relic neutralinos playing the role of a major dark matter 
constituent. 

The phenomenological analysis presented in this paper goes beyond the discussion of 
the experimental data specifically discussed here. We have tried to pin down the most 
relevant theoretical points, which are still at the origin of large uncertainties, and then 
require additional investigation. These are: i) size of the Higgs-quark-quark and the squark- 
quark-neutralino couplings, ii) properties of the WIMP distribution at low velocities (with 
the possible existence of a solar-bound WIMP population), iii) accurate determination of 
the propagation in the galactic halo and in the heliosphere for cosmic-ray antiprotons. 
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FIGURE CAPTIONS 



FIG. 1. Plot of ^ cricaiar'""^ veisus 171^. The solid line delimits the 3a C.L. annual- 
modulation region, obtained by the DAMA Nal(Tl) experiment with a total exposure of 
57 986 kg ■ day |0. This region was obtained by including the upper-bound constraints of 
Ref. 0, by setting pi at the standard reference value: pi = 0.3 GeV cm~^, and by taking 
into account uncertainties in the astrophysical velocities of the usual galactic Maxwellian 
distribution. Also shown in the present figure are the contour lines for the three values f o = 
170 km s~^ (short-dashed (red) line), Vq = 220 km s~^ (long-dash-short-dashed (blue) line), 
Vq = 270 km s~^ (long-dashed (green) line), separately. The scatter plot is calculated in the 
MSSM with the scan described in Sect. I; the points of the scatter plot are coded according 
to the value of the relic abundance, fl^h"^, of the relevant susy configuration: dots denote 
^x^"^ < 0.01, crosses denote 0.01 < fl^ < 0-1 empty circles denote > O-l- 

FIG. 2. Location of the DAMA annual-modulation region for four representative values 
of pf. pi = 0.1, 0.3, 0.5, 0.7 GeV cm~^. The scatter plots show only the configurations 
which lie inside the relevant annual-modulation region. A grey-level (color) code is used 
depending on the value of Vq employed in the extraction of the annual-modulation region: 
medium grey (red) denotes points which lie in the annual-modulation region extracted by 
setting vq = 170 Km s~^, dark grey (blue) denotes points which lie in the annual-modulation 
region extracted by setting vq = 220 Km s~^, light grey (green) denotes points which lie in 
the annual-modulation region extracted by setting vq = 270 Km s~^. The three sets are 
superimposed in that sequential order. 

FIG. 3. Scatter plot for the up-going muon flux from the center of the Earth for a 
standard Maxwellian distribution, Earth ^std^ versus m^. The grey-level (color) code is the 
same as in Fig. 2. The solid line denotes the 99.7% C.L. upper bounds, (<|)Earthyim^ derived 
from the data of the MACRO experiment 



FIG. 4. Enhancement effect in the up-going muon flux from the center of the Earth in 



case of a solar-bound population a la Damour-Krauss ^7j. The grey-level (color) code is 
the same as in Fig. 2. 

FIG. 5. Suppression effect in the up-going muon flux from the center of the Earth in 
case of the Gould-Alam conjecture [EH]. The grey-level (color) code is the same as in Fig. 
2. 



FIG. 6. Scatter plot for the up-going muon flux from the Sun for a standard Maxwellian 
distribution, versus m^. The grey-level (color) code is the same as in Fig. 2. The solid 

line denotes the 99.7% C.L. upper bounds, derived from the data of the MACRO experiment 
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FIG. 7. Scatter plot for the reduced Xr^ iii ^ comparison of the calculated cosmic-ray 



antiprotons fluxes with the combined experimental data of BESS95 and BESS97 ^J]. The 
horizontal line denotes the value Xr = 2.44, which for 13 d.o.f. corresponds to a 99.7% C.L., 
above which we disallow susy configurations. The grey-level (color) code is the same as in 
Fig. 2. 

FIG. 8. As in Fig. 2, once the constraints from the up-going muon fluxes from the center 
of the Earth are applied, assuming a Maxwellian halo distribution for relic neutralinos. The 
grey-level (color) code is the same as in Fig. 2. 

FIG. 9. Covering of the annual-modulation regions, if the constraint (<|)Earth^DK ^ 
^^Earthyim t^qj-q applied. The grey-level (color) code is the same as in Fig. 2. 

FIG. 10. Scatter plot for set T in the plane rrih - tan/?. The grey-level (color) code 
is the same as in Fig. 2. For each panel, the lower dashed line denotes the frontier of the 
complete scatter plot; the upper dashed line denotes the frontier, when only set 1 for the 
hadronic quantities of Sect. II is employed. The hatched region on the right is excluded 
by theory. The hatched region on the left is excluded by present data from LEP []TB[ and 



CDF |T^. The solid line represents the 95% C.L. bound reachable at LEP2, in case of non 
discovery of a neutral Higgs boson. 

FIG. 11. Scatter plot for set T in the plane mo - m^. The grey-level (color) code is 
the same as in Fig. 2. For each panel, the upper dashed line denotes the frontier of the 
complete scatter plot; the lower dashed line denotes the frontier, when only set 1 for the 
hadronic quantities of Sect. II is employed. 

FIG. 12. Neutralino relic abundance f^^^^ versus m^, once the constraints from up-going 
muon fluxes and cosmic-ray antiprotons are applied. The hatched region is disallowed by 
the upper limit on cold dark matter VLcDuh^ ^0-3 [iH- 
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FIG. 1. Plot of ^ cr^caiar^""'' veisus m^. The solid line delimits the 3a C.L. annual- 
modulation region, obtained by the DAMA Nal(Tl) experiment with a total exposure of 
57 986 kg ■ day |p. This region was obtained by including the upper-bound constraints of 
Ref. 0, by setting pi at the standard reference value: pi = 0.3 GeV cm~^, and by taking 
into account uncertainties in the astrophysical velocities of the usual galactic Maxwellian 
distribution. Also shown in the present figure are the contour lines for the three values vq = 
170 km s~^ (short-dashed (red) line), vo = 220 km s~^ (long-dash-short-dashed (blue) line), 
fo = 270 km s~^ (long-dashed (green) line), separately. The scatter plot is calculated in the 
MSSM with the scan described in Sect. I; the points of the scatter plot are coded according 
to the value of the relic abundance, Q^h"^, of the relevant susy configuration: dots denote 
flyh^ < 0.01, crosses denote 0.01 < fly h'^ < 0.1 and empty circles denote h? > 0.1. 



18 





50 100 15(3 200 EEO 100 150 aOO £B0 

(GeV) (GeV) 



FIG. 2. Location of the DAMA annual-modulation region for four representative values 
of pi: pi = 0.1, 0.3, 0.5, 0.7 GeV cm~'^. The scatter plots show only the configurations 
which lie inside the relevant annual-modulation region. A grey-level (color) code is used 
depending on the value of vq employed in the extraction of the annual-modulation region: 
medium grey (red) denotes points which lie in the annual-modulation region extracted by 
setting vo = 170 Km s~^, dark grey (blue) denotes points which lie in the annual-modulation 
region extracted by setting vq = 220 Km s~^, light grey (green) denotes points which lie in 
the annual-modulation region extracted by setting vq = 270 Km s~^. The three sets are 
superimposed in that sequential order. 
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FIG. 3. Scatter plot for the up-going muon flux from the center of the Earth for a 
standard Maxwellian distribution, (^<j)Earth-jstd^ versus m^. The grey-level (color) code is the 
same as in Fig. 2. The solid line denotes the 99.7% C.L. upper bounds, ($ Earth yim^ derived 



from the data of the MACRO experiment [45 
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FIG. 4. Enhancement effect in the up-going muon ffux from the center of the Earth in 



case of a solar-bound population a la Damour-Krauss The grey-level (color) code is 
the same as in Fig. 2. 
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FIG. 5. Suppression effect in the up-going muon ffux from tlie center of tlie Eartli in 



case of the Gould-Alam conjecture [^]. The grey-level (color) code is the same as in Fig. 
2. 
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FIG. 6. Scatter plot for the up-going muon flux from the Sun for a standard Maxwellian 
distribution, versus m^. The grey-level (color) code is the same as in Fig. 2. The solid 

line denotes the 99.7% C.L. upper bounds, derived from the data of the MACRO experiment 
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FIG. 7. Scatter plot for the reduced Xr^ iii a comparison of the calculated cosmic-ray 



antiprotons fluxes with the combined experimental data of BESS95 and BESS97 The 



horizontal line denotes the value Xr = '^■4:4:, which for 13 d.o.f. corresponds to a 99.7% C.L., 
above which we disallow susy configurations. The grey-level (color) code is the same as in 
Fig. 2. 
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FIG. 8. As in Fig. 2, once the constraints from the up-going muon fluxes from the center 
of the Earth are apphed, assuming a Maxwelhan halo distribution for relic neutralinos. The 
grey-level (color) code is the same as in Fig. 2. 
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FIG. 9. Covering of the annual-modulation regions, if the constraint (<|)Earth-jDK 
j^^Earthyim t^qj-q applied. The grey-level (color) code is the same as in Fig. 2. 
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FIG. 10. Scatter plot for set T in the plane rrih - tan/?. The grey-level (color) code 
is the same as in Fig. 2. For each panel, the lower dashed line denotes the frontier of the 
complete scatter plot; the upper dashed line denotes the frontier, when only set 1 for the 
hadronic quantities of Sect. II is employed. The hatched region on the right is excluded 
by theory. The hatched region on the left is excluded by present data from LEP |]TB[ and 
CDF The solid line represents the 95% CL. bound reachable at LEP2, in case of non 



discovery of a neutral Higgs boson. 
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FIG. 11. Scatter plot for set T in the plane mo - m^. The grey-level (color) code is 
the same as in Fig. 2. For each panel, the upper dashed line denotes the frontier of the 
complete scatter plot; the lower dashed line denotes the frontier, when only set 1 for the 
hadronic quantities of Sect. II is employed. 
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FIG. 12. Neutralino relic abundance f^^^^ versus m^, once the constraints from up-going 
muon fluxes and cosmic-ray antiprotons are applied. The hatched region is disallowed by 
the upper limit on cold dark matter flcDMh"^ ^0-3 [0- 
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